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Abstract 



We extract (for the first time) the ratio of the gluon condensate {g^ fatcG^) I {(^ sC^) expressed in terms of the liquid instanton radius 
Pc from charmonium moments sum rules by examining the effects of {OsG^} in the determinations of both and the running 
' ' MS mass mdmc). Using a global analysis of selected ratios of moments at different = 0, 4ml ™d 8m^ and keeping {a^G^) 
Ifrom 0.06 GeV"^, where the estimate of pc is almost independent of (a^G^), we deduce: pc - 0.98(21) GeV"' corresponding to 
'{g^fabcG^) - (31 + 13)GeV^(a5G^). The value of mdmc) is less affected (within the errors) by the variation of (ffjC^), where a 
CsJ [common solution from different moments are reached for {usG^) > 0.02 GeV^. Using the values of {a^G^} = 0.06(2) GeV"* from 
•some other channels and the previous value of {g^fabcG^), we deduce: mdmc) = 1261(18) MeV and mdmb) - 4173(10) MeV, 



where an estimate of the 4-loops iO(aj.)) contribution has been included. Our analysis indicates that the errors in the determinations 



-of the charm quark mass and of a, without taking into account the ones of the gluon condensates have been underestimated. To that 
CO 'accuracy, one can deduce the running light and heavy quark masses and their ratios evaluated at Mz, where it is remarkable to notice 
the approximate equalities: m^jmu » mh/m, ~ m,/mi, » 51(4), which might reveal some eventual underlying novel symmetry of 
the quark mass matrix in some Grand Unified Theories. 



.Keywords: QCD spectral sum rules, gluon condensates, heavy quark masses. 
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1. Introduction 

Non-zero values of the gluon condensates have been advocated 
by SVZ [1,2]. Indeed, the gluon condensates play an important 
role in gluodynamics (low-energy theorems,...) and in some 
bag models as they are directly related to the vacuum energy 
density (with standard notations): 

. /3{a,), 



Bar? 



■(ff,G2) 



(1) 



Moreover, the gluon condensates enter in the OPE of the 
hadronic correlators [1] and then are important in the analy- 
sis of QCD spectral sum rules (QSSR), especially, in the heavy 
quarks and in the pure Yang-Mills gluonia channels where the 
hght quark loops and quark condensates ' are absent to leading 
order [3-5]. The SVZ value: 



{a,G^} ^ 0.04 GeV"* 



(2) 



extracted (for the first time) from charmonium sum rules [1] 
has been challenged by different authors [3-5]. Though there 
are strong indications that the exact value of the gluon conden- 
sate is around this value or most likely 2 times this value as 
obtained from heavy quarks exponential moments [3-6], heavy 
quark mass-splittings [7] ande^e" [8-10], most recent determi- 
nations from T-decay [11-13] (see however [14]) and the pre- 
vious charmonium moments [15] indicate that its value is not 
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well determined. In fact, at present, the structure of the QCD 
vacuum is not yet under a good control. If one follows the SVZ 
idea based on the ordinary OPE, the QCD confinement can be 
parametrized by the sum of quark and gluon condensates of 
higher and higher dimensions ' . In order to estimate the higher 
dimension condensates, one usually assumes factorization us- 
ing vacuum saturation (leading 1/A^c approximation). However, 
in many examples, this assumption appears to be badly vio- 
lated [8, 9, 11-13, 18-24]. Different phenomenological works 
have been performed for understanding the complex structure 
of the QCD vacuum in the V + A and V - A channels of the 
light flavours [8, 9, 1 1-13, 18-24] and from lattice calculations 
[25, 27, 28]. Here, we shall estimate (for the first time) the ra- 
tio of the dimension-6 {g^fatcG^) over the dimension-4 {asG^} 
gluon condensates using charmonium sum rules ^, in the aim 
to clarify the different inaccurate proposals from some instan- 
ton liquid models. In so doing, we find that it is convenient to 
introduce the instanton radius pc'. 



g'fahcG'} 



4 12;r 
5~pf 



(3) 



The value of pc ranges from 1/3 fm=1.5 GeV ' [29], 0.5 
fm=2.5 GeV-i [30] to 0.9 fm= 4.5 GeV"' [1]. As (g^fabcG^) 
contributes like 1/p^ in the OPE analysis, a more precise value 
of Pc is crucial for checking the convergence of the OPE. The 



' A possible existence of an additional l/Q^ term induced by large order 
terms of PT series has been discussed in [16, 17]. 

^The (g^fabcG^) condensate does not contiibute in the chiral limit = 
in the vector and axial- vector channels of Ught flavours. 
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estimate of pc from charmonium sum rules is feasible as the 
light quark condensates niqiqq) contributes, to higher loop or- 
der and which are chiral suppressed are negligible, while the 
heavy quark condensate contribution can be absorbed into the 
gluon one as mentioned earlier. 

2. Moment sum rules 

Here, we shall be concerned with the two-point correlator of a 
heavy quark Q: 



the Appendix of [15] to 3-loops ((3(a;^)) accuracy in terms of 
the running heavy quark mass ^ using the pQCD results of [32], 
while the 2^ = moments to 4-loops {O(a^J) are given in [33] 
using the pQCD results in [34, 35] Among the different mo- 
ments given in [15], we shall select three moments where both 
the (a.,)" (n - I, 2), the gluon condensate contributions and 
the effects of the higher resonances plus the QCD continuum 
are relatively small but not negligible. These conditions can be 
simultaneously satisfied by the moments ^: 



At2,3,4 
At8,9,10 



for 
for 



2=0, 
= Ami 



d'x e-^'^^OlTJ^Qix) (4(0)) 10) , (4) j^^^^^^^^^ fo^ = Sm^ . 



(8) 



where : Jq - Qy^Q is the heavy quark neutral vector current. 
Im IlcCi) can be related to the charmonium leptonic widths and 
masses. In a narrow width approximation (NWA): 



3; and Y^- 



(5) 



where Nc = 3; and Yip^e+e- are the mass and leptonic width 
of the y/i/r mesons; Qc = 2/3 is the charm electric charge in 
units of e; a = 1/133.6 is the running electromagnetic coupling 
evaluated at M^. We shall use the experimental values of the 
J/i// parameters compiled in Table 1. 

Table 1: Masses and electronic widths of the J/ili family from PDG 08 [31], 



Name 


Mass [MeV] 


[keV] 


jms) 


3096.916(11) 


5.55(14) 


i/r(2S) 


3686.093(34) 


2.33(7) 


iA(3770) 


3775.2(1.7) 


0.259(16) 


iA(4040) 


4039(1) 


0.86(7) 


i/'(4160) 


4153(3) 


0.83(7) 


1^(4415) 


4421(4) 


0.58(7) 



One may also work with more moments but these will not bring 
newer informations. Lower moments are more sensitive to the 
experimental errors and to the QCD continuum while higher 
moments are more sensitive to higher dimension condensates 
which are not under a good control''. Moreover, one can also 
note from the QCD expressions given by [15] that for n larger 
than in our previous selected choice, the signs of the pQCD 
corrections start to change compared to the original ones of the 
two-point correlator. A such change may introduce some sys- 
tematical difficulties inherent to the approach for controlling the 
size of higher order terms ^. 
We shall work with the ratios of moments: 



r2/3 and 
rg/g and 
ri3/i4 and 



1-2/4 for 
n/io for 
n3/i5 for 



2=0 

2' = Ami 
= Sml 



and use as inputs, in this first step: 
mdmc) = 1.26(3) GeV, 



(9) 



(10) 



as given by different approaches using charmonium moments 
sum rules [1, 3-5, 15, 31, 33, 41-43] and which we shall re- 
estimate later on. We shall also use: 



Different forms of QSSR exist in the literature [3-5]. We shall 
work here with the moments: 



^ ' J 4ml *■ 

and more likely with their ratios: 
M„ 



Kit, ml) 



r„/„+iiQ ) = 



Al„+i' 



r„/n+2iQ) = 



Mn_ 

M„+2 ' 



(6) 



(7) 



where the experimental sides are more precise than the abso- 
lute moments M„. Also, in the ratios, partial cancellations of 
different perturbative as well as non-perturbative terms occur, 
which render the QCD approximation more precise than in the 
absolute moments. 

The QCD sides of the sum rules are known in the literature 
since the original works of SVZ [1]. Their expressions at the 
subtraction scale = are given explicitly numerically in 



a.iMr) = 0.3249(80) =^ Qr,(mc)L^=4 = 0.408(14) (11) 



from T-decay [14]; a value which agrees with the central value 
of the world average [3 1, 44] when rurmed until Mz. 
The QCD expressions of the moments and their ratios are given 
in Table 2. 



'We shall use these expressions in our analysis and we shall correct our 
final results on the quark masses by adding an estimate of the 4-loops {0{al)) 
contributions. 

'*Some Pade approximants are given in [36]. 

'However, one should note that the accuracy of the = moments is less 
than that of the ?i moments, while one cannot use higher moments due to 
the bad convergence of the OPE in this case. 

^The contributions of the dimension-8 condensates have been evaluated in 
[37] and can be sizeable if one assumes factorization which might not be ap- 
phed here [38]. 

'in [39], some low energy gluon contributions to order o-^ to the correlator 
can invahdate the uses of the = 0-moments for n > 4 (see however [40]). 
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Table 2: QCD expressions of the moments and their ratios normalized to the lowest 
order terms and given to order a; and including the dimension-4 and -6 gluon condensates 
deiived from [15]; d4 = (a,,G^)/(4mJ)^ ^ (1.49 ± 0.50) X 10"' if we use = 1.261 
GeV from Table 3 and <Q-.,C^> = 0.06 GeV" in Eq. 17; p64 = (g- fcbcC^)H(a,G'^)Aml) =a 
(4.88 ± 2.05) is the ratio between the dimension-6 and -4 gluon condensate contributions 
if we use the numerical value in Eq. 13; = a jn. 



Mom 



QCD expression 



Ml 

Mi 

M4 

ri/i 
1/2 
'2/4 

Ms 

M9 

Mio 

rs/9 
1/2 
''s/io 

Q' = - 

Ml4 

Mis 

1/2 

r 

13/15 



1 + 2.427a, + 6.110a? - ^4 (18.61 - 0.83p64) 



1 + 1.917a, +6. 115a 
1 + 1.100a, + 4.402a 
1 + 0.510a,, -0.983a 
1 + 0.664a, - 0.095a 



- d4 (45.71 -4.00p64) 

- d4 (90.21 - 12.76p64) 
+ d4(27.10-3.17p64) 
+ rf4 (35.80 - 5.97p64) 



4ml 



1 + 1.118a, + 4.253a2 - (77.51 - 5.02p64) 
1 + 0.601a, + 2.700a2 - (104.86 - 9.20p64) 
1 + 0.045a, + 1.136a2 - d^ (137.83 - 15.63p64) 
1 + 0.517as + 1.242a2 + ^4 (27.35 - 4.18p64) 
1 + 0.537a,, + 1.390aJ + d^ (30.16 - 5.31p64) 



1 + 0.776a, + 3.061a2 - ^4 (90.37 - 5.19p64) 



1 +0.4 12a 
1 +0.03 la 
1 + 0.364a, 
1 + 0.373a 



+ 1.909a2 - ^4 (109.75 - 7.89p64) 
+ 0.770a2 - d^ (137.72 - 1 1.53p64) 
+ 1.002a2 + d4 (19.38 - 2.70p64) 
+ 1.065a2 + d4 (23.68 - 3.17p64) 



3. pc from charmonium ratios of moments 

We shall work with the ratios of moments in Eq. 9. We 
parametrize the spectral function by a sum of the six 7/i^-like 
narrow resonances below 4.6 GeV and use its pQCD expres- 
sion from yfi = (4.6 ±0.1) GeV. We extract the value of for 
a large range of {asG^}. 

One can see in Fig. 1 that the results are very stable for all 
moments for (a^G^) > 0.06 GeV"*, from which we deduce the 
value of Pc given in Table 3. The errors in pc come respec- 
tively from the values of nic and of the choice of the moments 
at given Q~. These errors are included in the regions given in 
Fig. 1. The one due to a, and to experiments are negligible . 
Our final averaged result is: 



p, = 0.98(21) GeV" 



for (a,G^) > 0.06 GeV"* 



(12) 



which we consider as an improvement of the different estimates 
based on instanton liquid models [1, 29, 30] recalled in the in- 
troduction. However, our result agrees within the error with the 
one in [29] but is smaller by a factor of about 4 than the SVZ 
estimate [1]. Using Eq. 3, our result for pc corresponds to ^: 

{g'fabcG'} = (31 + 13) GeV\a,G^} , (13) 



'^One can improve this parametrization by taking into account finite width 
corrections using BES data [45], but these corrections will be negligible in the 
moments which we shall use. 

'a recent estimate using exponential sum rules and including the 
dimension-eight condensate leads to a lower value (g^fabcG^) = (8.3 ± 
1.0) GeV2<ff,,G-> [46], which is still higher than the previous dilute gas in- 
stanton estimates. 
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Figure 1 : Determinations of pc in GeV ^ versus {a^G'^) in GeV^ from different 
sets of moments : Q- = (red: dashed-dotted) ; = 4m^ (green: continuous); 
= Hm^ (blue: dotted). 



Table 3: Ratio of the {g^ fahcG^)l{<^sG~) as defined in Eq. 3, and value of m,-(mc) 
from charmonium moments known to 3-loops. The errors in p^. come from the choice of 
moments and from the en"or on given in Eq. 10. The ones due to tr^, {a^G-) and to the 
data on the Jjip family are negligible. The value of nic is taken at {OsG^) - 0.06(2) GeV* 
as given in Eq. 17. The eiTors on come respectively from the ones of pc and (or^G"). 
The ones due to tr^ and the data are about 1 MeV each which are negligible. 



Mom 


Pc [GeV-'l 


m^imc) [MeV] 


e-=o: 

''2/3 > ''2/4 


0.800(490) 


1234(34)(8) 


2-=4mJ : 

''8/9^ ''8/10 


1.025(425) 


1265(29)(9) 


e-=8m2 . 

''13/14. ''13/15 


1.025(275) 


1268(18)(7) 


Average 


0.98(21) 


1261(15) 



indicating that it is much bigger than usually assumed in the 
literature. It is also smaller than the lattice result in S U(2) pure 
Yang-Mills [25] and than a rough estimate extended to S U (3) 
with dynamical fermions [26] using the result in [25, 27]. One 
should notice that this value of the gluon condensate has been 
extracted by assuming that the OPE including the dimension-6 
gluon condensate gives a good description of the experimental 
data which implicitly assumes that the contributions of higher 
dimension condensates are negligible in the analysis. In other 
words, one may also interpret this value as the one of an "effec- 
tive gluon condensate " which parametrizes all higher dimen- 
sion condensates contributing to the OPE. Fitted values of the 
higher dimensions vacuum condensates in the light quark chan- 
nels have been also found to be larger than the vacuum satu- 
ration assumptions [8, 9, 18-20, 24] and our results seem to 
go towards this direction. However, as we shall expUcitly dis- 
cuss later, these values of the gluon condensates remain still a 
correction compared to the one of the lowest order perturbative 
contribution in the ratios of moments which we use here and do 
not break the OPE. 
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4. mdmc) from charmonium ratios of moments 
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Figure 2: Determinations of mdnic) in GeV versus {a,G'^) in GeV* from dif- 
ferent sets of moments . = Q (red: daslied-dotted) ; (fi = 4ml (green: 
continuous); Q- = 8m^ (blue: dotted). 



For extracting mdmc), we equate the QCD and experimental 
sides of the ratios of moments r„i„+\{Q^) and r„/„+2 (g^) and 
solve exactly numerically these equations using the Mathemat- 
ica subroutine FindRoot. Though the equations look simple for 
the - moment: 



c d e 
F,h(x) = ax + b + - + ^ + —r ^ F,,,p(x) 



(14) 



X = mf., they are highly non-trivial for - 4-m^, 8m^ and for 
higher values of n due to the appearance of the term: 



8,/, 



(Ml + Q^r 



(15) 



in the experimental side of the ratio of moments. We show in 
Fig. 2 the different solutions of mdmc) versus {osG^) for each 
ratio r„/„+i((2^) and rnin+dQ^) of Q^-moments. One can notice 
that common solutions of different ratios moments occur in the 
range of values: 



{a,G^) > 0.02 GeV* , 



(16) 



indicating that the central values of mdmc) are not very sensi- 
tive to the one of {osG^}. This feature confirms the unconclu- 
sive range of values obtained for {UsG^} in [15]. However, zero 
and negative values of {a^G^) as obtained from some analysis 
of T-decays [13] '"are excluded by our present result and by the 
one in [15]. The slight difference between [15] with our analy- 
sis is that we put implicitly the previously determined value of 
{g^fabcG^} and its correlation with {UsG^) in the extraction of 
mdmc). This fact explains the increase of the errors for increas- 
ing values of {asG^) in our analysis. In order to improve the 
determination of mc, we use the values: 



{a,G^) = 6(2) X IQ-^ GeV"* 



(17) 



'"One sliould notice tliat due to the kinematical structure of tlie original r- 
decay widtli [47], tlie gluon condensate contribution acquiers there an extra ff , 
coefficient compared with the one of the two-point correlator which suppresses 
its contribution and can render inaccurate its extraction from ffiis observable. 



obtained by enlarging the error of the average value 0.06(1) 
GeV"* from the heavy quarkonia mass-splittings [7]: 



{OsG^) = 7.5(2.5) X IQ-^ GeV* , 
and from e^e" — > / = 1 hadrons sum rules [10]: 
{asG^) = 6.1(0.7) X IQ-^ GeV* . 



(18) 



(19) 



These previous values agree with the one about 0.069 GeV* ob- 
tained from S U(3) lattice with dynamical fermions [27]. Using 
this value, we show in Table 3 the value of mdmc) from each 
sets of ratios of moments from which we can deduce the mean 
value: 



mdmc)\3-ioops ^ 1261(15) MeV, 



(20) 



obtained from a 3-loop {0(a^^)) expression of the ratios of 
moments. This value is more weighted by the one from the 
- 8m^ ratios of moments which give the most precise pre- 
dictions. 



5. Comments on the results 

•2^ — moments 
We note that our result mdmc) = 1234(35) MeV from the 
- moments agrees within the error with the four-loops 
recent estimate 1279(13) MeV in [33] based on the lowest 
n - I moment. Our result is less precise due mainly to 
the errors induced by the presence of {g^fabcG^) and of its 
correlated {usG^) condensate in our analysis and to the error 
induced by the choice of the ratio of moments as can be seen 
in Fig. 2. It is informative to compare the size of each QCD 
corrections in the OPE, that can be deduced from the QCD 
expressions given in Table 2: 

- One should first notice that each perturbative and non- 
perturbative corrections tends to partially cancel out in the 
ratios of moments, which render the QCD PT series and OPE 
more convergent for the ratios than for the corresponding 
individual moments. Therefore, one expects that these ratios of 
moments can lead to more robust predictions. 

- Using the previous values of the QCD parameters and for 
definiteness {a,G^) = 0.06 GeV*, one finds that the contribu- 
tion of {osG^) in MiiQ) is about -2.8% which is comparable 
with the one -1.6% from -Q.lZa] [33], while, for Md^), it 
is -6.8% which is about the one 10% from a]. For Al2(0) the 
ig^fahcG^} contribution is about 0.6% which is about 1/3 of the 
a] one, while for A43(0) it is about 3% compared with 2% from 
0.299aj and with 10% from a^. These features indicate that 
the non-perturbative corrections can be comparable with the 
PT radiative corrections and cannot be neglected like usually 
done in the exisiting Uterature (see e.g. [33] and references 
therein). The same remark also applies to the extraction of as 
in [48] from low-« moments. 

- Finally, the leading experimental error due to the J/t// leptonic 
widths, which gives a strong limitation to the accuracy of the 
low n moments (see e.g. [33] and references therein), partially 
cancel out in the ratio of moments such that the experimental 
error induces only a negligible error of about 1 MeV in the 
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dermination of nic. In the same way, the high mass states 
contributions to the spectral functions are more suppressed 
in the ratios of moments, which then avoid some difficulties 
induced by the present data in the high-energy regions. 

• g2 Q Yfioments 

With these moments, we can work at larger values of n, where 
the experimental sides of the sum rules become more accurate 
due to the increase of the weight of the lower mass resonances 
contributions in these moments: 

- Comparing e.g. the QCD sides of the n = 2, = MiiG) 
with that of the n - - Am^ A^8(4m^)-moments, we find 
from Table 2 that the sum of the PT corrections up to order 

for Ms(4ml) (21%) are about 1/2 of the ones for MiiO) 
(42%). Moreover, though the size of the sum of the NP terms 
increases from -2.2% for Al2(0) to -9.6% for M^{Aml), one 
can see that the ratio between the (g^fabcG^) over the {a,G^) 
contributions are almost unchanged of about 20%, indicating 
the good convergence of the OPE in the analysis. 

- One can also deduce from Table 2, that the PT and NP QCD 
corrections are much lower for the ratios of moments. The sum 
of PT corrections is typically 6.5% while the NP ones are 2%. 
A convergence of the OPE is still observed though the ratio 
between the {g^fatcG^} over the {(XsG^) contributions is larger 
(0.6 ~ 0.8) for the ratios of moments than for the moments 
(0.3 ~ 0.4). 

- Our results from these + ratios of moments given in 
Table 3 agree with the ones 1275(15) MeV obtained in [15] 
though less accurate due to the effect of the error on (g^fabcG^) 
included in our analysis. 

• Concluding remarks 

- These previous facts indicate that a precise determination of 
mdnic) and of or^ requires the inclusion of the non perturbative 
condensates which can induce large errors even for the lowest 

=Q moments and which have not been taken properly into 
account in the existing Uterature. 

- One can also note from Table 3 that the best estimate of 
nicimc) comes from the ratio of moments where the 
error due to the choice of the ratio of moments is smaller than 
in the case =0. 



which we consider to be more conservative than the one of 
about 2-3 MeV given in [15]. Alternatively, one can also 
minimize the v-dependence by working at large and with 
low n moments and after running down the result to nic. In this 
way, one would obtain a slightly smaller error of about 3-5 
MeV [33]. 

• Shift due to Coulombic corrections 

The contribution due to Coulombic corrections are expected to 
be negligible (about 1-2 MeV [15]) because the system is still 
relativistic. In fact, the Coulomb radius: 



^Coul 



mcCpasimc) 



^ 3 GeV 



(23) 



{Cf - 4/3) is much larger than the confinement radius rconj ~ 
1 GeV ' . These corrections can even be made much smaller by 
working with a 2^ moments rather than with aQ^ = one, 
as the quark velocity behaves for large n as: 



+ QV4m^) /n , 



(24) 



which, e.g., for = 8m^ and n = 14, is about 0.46. This value 
is not small and not inside the nonrelativistic region. We can ap- 
proximately estimate this effect by working with the resummed 
Coulombic expression of the spectral function [49]: 



Coul 



2 1 



(25) 



where: x = nCpas/v, Cp = 4/3 and v = ^1 - 4mg/f. We com- 
pare the value of the ratio of moments using this perturbative 
expression for the spectral function with the one obtained from 
PT theory including radiative corrections. In the case = 8m^ 
and n = 14, where the most precise result is obtained, the cor- 
rections induced by the Coulombic contributions to the value of 
nic are negligible " : 



^mJc««/ = -(0.4±0.4)MeV, 



(26) 



where we have assumed that our determination is known within 
100% error. 



• Error due to the subtraction point v 

Our previous results in Eq. 20 have been obtained at the 
subtraction point = w^. If the whole series is known, the 
results should be independent of v. The knowledge of the PT 
series to 4-loops decreases the sensitivity of the results on v. 
To order where the moments have been evaluated, one can 
introduce this v-dependence through the replacement (see e.g.: 
[3, 5]): 

a.imc) ^ a,(v)x(l-A ^log— ) , (21) 

where P\ - -(1/2)(11 - 2nf/3) for n/-flavours. Taking 0.5 < 
v^/ml < 2 , we deduce from the = 8m^ ratio of moments: 



6. mdmc) to order Oio^) 

• Estimate of the 0(al) and higher order corrections 
Observing that the coefficients of the PT corrections for the mo- 
ment decrease when n increases and do not flip sign compared 
with the lowest moments (see Table 2) and assuming that the 
ratio of the 3-loop over the 4-loop coefficients are approxima- 
tively the same for each moments, we can write to 4-loops: 

Mi3(8m2) ~ 1 + 0.78fl, + 3.06fl? - 5.6fl^ 
Mui^ml) ~ 1 -H 0.41a, -H 1.91a^-3.5a^ 
Mi5(8OTc) ~ 1 + 0.03a, -I- 0.77a? - 1.4a^ , (27) 



±6 MeV , 



(22) 



' ' Some further arguments justifying the smallness of these contributions can 
be found in [15]. 



5 



where a, = ajn. We have used the 4-loops coefficient -5.6 
obtained in [34, 35] for low n = I, = moment, which 
we expect to be an overestimate of the coefficient of A1i3(8m^). 
These expressions lead to an corrections of about 2% for 
both ri3/i4(8m^) and ri^/isiSm^) and an correction of about 
-0.3% and -0.9%. Then, taking the average of the two correc- 
tions, we may expect that the a] corrections can provide a max- 
imal shift of the charm quark mass of about -0.3% leading to 

Sm, U-ioops ^±(2x4) MeV , (28) 

a range of values expected from some alternative estimates 
[33]. The factor 2 in front assumes ffie estimate of higher 
order PT (0(a") : n > 4) contributions or by duaUty ffie 1/s 
corrections due to ffie tachyonic gluon mass [16]. 

• Final value of Trie to order Oia^^) 

Adding ffie previous estimates of new sources of contributions 
and errors into ffie 3-loops result in Eq. 20, we obtain to 0{al) 

mc(mc)\4-ioops = 1261(18) MeV . (29) 

This result is comparable with ffie existing ones obtained from 
moment sum rules in the literature [1, 3-5, 15, 31, 33, 41^3]. 
Our final result confirms and improves (reduction of errors) ear- 
her sum rules analysis obtained using PT lower orders charmo- 
nium sum rules [1, 41, 42] It is in agreement wiffi ffie most 
recent results from sum rules [15, 33] mentioned previously 
with the lattice determination 1268(9) MeV [50] and wiffi ffie 
PDG 08 average [31] (l.27+[;;'^]) GeV. 

7. Determination of mgCm^) 

We extend the previous analysis of the charmonium system to 
bottomium. In the following, we shall use ffie value: 

aAmb)\nf=5 = 0.219(4) , (30) 

deduced from asirrir) in Eq. 1 1 . We shall use as experimental 



Table 4: Masses and electronic widths of the T family from PDG 08[31]. 



Name 


Mass [MeV] 


Tt^e*e- [keV] 


Y(15) 


9460.30(26) 


1.340(18) 


Y(2S) 


10023.26(31) 


0.612(11) 


Y(35) 


10355.2(5) 


0.443(8) 


Y(4S) 


10579.4(1.2) 


0.272(29) 


Y(10860) 


10865(8) 


0.31(7) 


Y(11020) 


11019(8) 


0.13(3) 



More complete references can be found in Table 53.5 page 602 of [3]. 

However, one should mention that a sum rule analysis of the D meson 
mass using the pseudoscalar correlator to order leads to [52] 'tnc(mc)\pseudo = 
1 . 1 0(4) GeV but the value of fo, agrees with the present lattice calculations. We 
shall reconsider this point elsewhere. 



inputs ffie T-family parameters in Table 4 using NWA and 
parametrize ffie spectral function above V? = (11.098 ± 0.079) 
GeV by its pQCD expression (QCD continuum), where the 
error in the continuum ffireshold is given by the total widffi 
of ffie T( 11020). Using ffie previous moments, ffie dominant 
contributions will come from the two lowest ground states 
while finite widffi corrections will not be observable. We show 
in Table 5 the results from different moments known to 3-loops. 



Table 5: Value of mi,(mh) from bottomiun moments known to 3-loops. The errors on 
m/, come respectively from the choice of the moments, a^, the data on the T family and 
the choice of the QCD continuum threshold. The ones due to the gluon condensates are 



negligil^le here. 


Mom 


nibimb) [MeV] 


2^=0: 




''2/3! ''2/4 


4160(4)(2)(3)(3) 


e'=4m^ : 




?'8/9, ''8/10 


4177(2)(3)(3)(6) 


Q^=?.ml : 






4183(2)(4)(2)(6) 


Average 


4173(4) 



• Error due to the subtraction point 

We study the effect of the subtraction point by taking 
0.5 < v^/m^ < 2 and using the expression in Eq. 21. We 
induce an error: 

5^,1, = ±6 MeV. (31) 

This error can be furffier reduced by using ffie 4-loops expres- 
sion of the moments. 

• Shift due to Coulombic corrections 

Using nihimb) - 4.24 GeV into Eq. 23, one obtains: 

rco„;- 1.6 GeV-^ (32) 

which is still larger than ffie confinement radius rconf ~ 1 
GeV"^ . These corrections can be render much smaller by work- 
ing with a 2^ = 4m^ and large n moments rather than with a 
- one, where ffie fe-quark velocity is about 0.45 from Eq. 
24. This value is still inside the relativistic region, where one 
can safely neglect these Coulombic corrections. Indeed, using 
ffie previous expression of ffie Coulombic corrections in Eq. 23, 
we obtain the shift: 

6m,\coui ^ ±6MeV , (33) 

which is about ffie same value as ffie one obtained in [43]. 
Due to the theoretical uncertainties on the real effect of the 
Coulombic corrections, we consider our previous estimate as 
another source of errors but not as a safe correction. 

• Value ofnibimb) to 0{al) 
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In the case of the b quark, our previous estimate of the 4-loops 
{0{a^)) contribution induces an error: 



Sm„U-i„op. - +(2x2) MeV, 



(34) 



where the factor 2 is assumed to include higher order or/and 
/l^-tachyonic gluon mass corrections. Adding these new cor- 
rections to the one in Table 5, we deduce to 4-loops accuracy: 



mt,{mb)\A-ioops = 4173(10) MeV , 



(35) 



which is relatively more precise than that of nic as the non per- 
turbative contributions are much smaller here, while is eval- 
uated at a higher scale. This result is in excellent agreement 
with the one: 



nihimb) = 4171(14) MeV, 



(36) 



obtained in [51] using ratios of moments based on the criteria 
of stabilities versus the degree n (number of 2^-derivatives) of 
moments and including a] and the dimension eight condensates 
contributions. This result is also in good agreement with with 
the PDG average [31] : 



mMpDG = (4.20!°:J^) GeV , 



(37) 



but more precise. It also agrees with some previous results 
quoted in Table 53.6 (page 603) of [3]. However, it is worth 
mentionning that hke in the case of the Z)-meson, the analysis 
of the B meson mass from the pseudoscalar sum rule to order 
leads to a lower value of (4.05 ± 0.06) GeV [52] which will 
be reconsidered elsewhere. 



8. Running light and Iieavy quark masses at Mz 

For direct uses in some phenomenological applications and as 
inputs in some Grand Unified Model Buildings, it can be use- 
ful to convert these running masses mgimQ) to the ones eval- 
uated at the Z-mass. This can be easily done after taking care 
on different quark threshold eff'ects. Using, e.g., the Mathe- 
matica RunDec package [53], we deduce, to 4-loops accuracy, 
from Eqs. 29 and 35, the rurming masses evaluated at Mz for 5 
flavours: 



m,(Mz) = 616(9)^„(6)„, MeV , 
mMz) = 2920(7)„,(22)„, MeV . 



(38) 



The errors are due respectively to the values of the running mass 
and of Qs when one performs the QCD evolutions. 
In a similar way, we can also deduce the ones of rhu4,s{Mz) 
and m,(Mz) by using respectively the average value from QSSR 
predictions to 4-loops[54, 55]: 



m,(2) = 96.1(4.8) MeV, 
mil) = 5.1(2) MeV, 
m„(2) = 2.8(2) MeV , 



(39) 



and the on-shell top quark mass average to 3-loops [31]: 

= 171.2(2.1) GeV. (40) 



We obtain for 5 flavours: 

m,(Mz) = 53.9(2.9)„X1.9)a,MeV. 
rnaiMz) = 2.47(10)^/3)„, MeV , 
m„(Mz) = 1.30(10)„„(3)„.MeV, 



and: 



m,iMz) = 168.4(2.1)„,(0.1)„. GeV . 



(41) 



(42) 



Combining the previous results, we obtain the ratios of rurming 
masses at Mz: 



m 

rric 



■ 4.7(1) , 



(43) 



and: 



: 42(5) , 



58(2) , 



(44) 



Mh m, 
— = 54(4) , — 

OT„ Ms nib 

where it is remarkable to notice that the ratios nis/mu, Mb/ Ms 
and mi/mi, are almost equal which might reveal some eventual 
underlying novel symmetry of the quark mass matrix [56, 57] 
in some Grand Unified Theories [58]. One should also observe 
that when one runs the ratio ms/nig from 2 GeV to Mz, the 
central value is not strictly constant (contrary to what expected 
from its renormalization group invariance) though the two val- 
ues agree within the errors. This is due to different threshold 
effects and to the truncation of the series at a given order of PT. 

9. Conclusions 

- Firstly, our analysis has been motivated to extract (for the first 
time) from the sum rules, the instanton liquid model radius pc 
in Eq. 12 , which parametrizes the ratio of the (g^fabcG^) over 
the {a^G^) gluon condensates, where the corresponding value 
of ig^fabcG^) in Eq. 13 is much larger than usually quoted in 
the literature. However, despite this large value, the OPE in 
the ratios of moments which we have used continues to present 
a good convergence. Because we have neglected the contri- 
butions of higher dimension condensates in our approach, one 
may consider this value of {g^fahcG^} as that of an "effective 
gluon condensate" which may include in it all the higher di- 
mension condensates contributing to the OPE and not consid- 
ered in our analysis. 

- Using the previous result and thanks to the recent progresses 
in evaluating accurately the pQCD series of the heavy quarks 
vector correlators [32, 34, 35] and to more accurate measure- 
ments of the corresponding spectral functions [31, 45, 59], it be- 
comes possible to extract with a high precision the heavy quark 
masses using higher n ratios of moment sum rules. The results 
for nic and nib in Eqs. 29 and 35, where different sources of er- 
rors are under a good control (see also the comments in section 
5), are among the most accurate measurements available today. 
These results confirm and improve estimates done in the early 
days of sum rules [1, 3-5, 41, 42] 



'^One can notice, in different papers written by the author for extracting ntc 
and nib from heavy quarkonia sum rules, that their central values remain very 
stable since the 1st paper in 1987 [42]. This feature indicates (a posteriori) the 
self-consistency of the approach and the good convergence of the PT and OPE, 
especially when one works with the running MS masses. 
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- Compared with some other recent determinations based on 
low-n and =0 moments [32-35] our approach should (a 
priori) be more accurate because we work with ratios of mo- 
ments which are less sensitive to the continuum contribution 
than the individual moments. We also refrain to take too high 
moments where their QCD expressions can become difficult to 
control. The apparent accuracy of the results obtained in the 
current literature are also due (among others) to the neglect of 
the gluon condensates (see Fig. 2) which are one of the main 
sources of the errors in the determinations of rric (and of as) 
from moment sum rules (see Table 3). The same remarks also 
apply to the higher-n moments used in [15, 43]. 

- The agreement of the present results with the most precise 
recent lattice calculation of ntc [50] confirms the robustness of 
the higher n ratios of moments sum rule approach and, in gen- 
eral, the ability of QSSR to extract reliably the QCD parameters 
from hadron properties. It becomes now challenging to check 
our estimate of mt given in Eq. 35 using lattice calculations. 

- Finally, the approximate equalities of the diff'erent ratios of 
the quark masses in Eq. 44, when they are evaluated at Mz with 
the same number of flavours might reveal some eventual under- 
lying novel symmetry of the quark mass matrix in some Grand 
Unified Theories. 
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